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This  paper  summarises  a  literature  survey  on  high  temperature  measurement  and 
describes  the  design  of  probes  used  in  plume  measurements.  There  were  no  cases 
reported  of  measurements  in  extreme  environments  such  as  exist  in  solid  rocket 
exhausts,  but  there  were  a  number  of  thermocouple  designs  which  had  been  used  under 
less  extreme  conditions  and  which  could  be  further  developed. 


Tungsten-rhenium{W-Rh)  thermocouples  had  the  combined  properties  of  strength  at 
high  temperatures,  high  thermoelectric  emf,  and  resistance  to  chemical  attack.  A 
shielded  probe  was  required,  both  to  protect  the  thermcKouple  junction,  and  to 
minimise  radiative  heat  los.ses.  After  some  experimentation,  a  twin  shielded  design 
made  from  molybdenum  gave  acceptable  results.  Corrections  for  thermal  conduction 
losses  were  made  based  on  a  method  obtained  from  the  literature.  Radiation  losses  were 
minimised  with  this  probe  design,  and  corrections  for  these  losses  were  too  complex  and 
unreliable  to  be  included. 
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Design  of  Thermocouple  Probes  for 
Measurement  of  Rocket  Exhaust  Plume 
Temperatures 


EXECUTIVE  SUMMARY 

There  is  consider^hle  difficulty  and  expense  involved  in  tlie  n:\easurement  of  tiie 
infra-red  signature  of  rocket  and  aircraft  exhaust  phunes  under  in-service 
conditions,  and  it  is  often  impossible  to  measure  the  plumes  of  enemy  planes  or 
missiles.  The  alternative  is  to  calculate  the  appropriate  emissions  using  sophisticated 
computer  codes.  Measurement  of  temperatiue  profiles  of  rocket  exhausts  offers  an 
important  method  for  validating  these  codes. 

Temperature  profiles  can  be  measured  using  non-intrusive  methods,  but  these  are 
complex  and  expensive.  Intrusive  probes  based  on  thermocouples  are  a  possible 
alternative,  but  a  rocket  exhaust  presents  a  very  extreme  environment  of  high 
temperature,  high  vel<xity,  high  turbulence  and  chemical  reactivity.  These  factors 
place  severe  constraints  on  the  type  and  design  of  probe  that  can  be  used. 

This  paper  summarises  a  literature  siuvey  on  high  temperature  measurement  and 
descril^  the  design  of  probes  used  in  plume  measurements.  There  were  no  cases 
reported  of  measurements  in  extreme  environments  such  as  exist  in  solid  rocket 
eT^usts,  but  there  were  a  number  of  thermocouple  designs  which  had  been  used 
imder  less  extreme  conditions  and  which  could  be  fiurther  developed.  The  literature 
survey  also  covered  discussions  of  corrections  for  conduction  and  radiation  losses. 

Timgsten-rheniiun(W-Rh)  thermocouples  were  chosen  because  they  had  the 
combined  properties  of  strength  at  high  temperatures,  high  thermoelectric  emf,  and 
resistance  to  chemical  attack.  A  shielded  probe  was  used,  both  to  protect  the 
thermocouple  junction,  and  to  minimise  radiative  heat  losses.  After  some 
experimentation,  a  twin  shielded  design  made  from  molybdemun  gave  acceptable 
results.  Corrections  for  thermal  conduction  losses  were  made  based  on  a  method 
obtained  from  the  literature.  Radiation  losses  were  minimised  with  this  probe 
design,  and  corrections  for  these  losses  were  too  complex  and  unreliable  to  be 
included. 
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1.  Introduction 


It  is  very  difficult  and  expensive  to  measure  the  IR  sigiiature  of  rocket  and 
aircraft  exhaust  plumes  under  service  conditions,  and  it  is  often  impossible  to 
measure  the  plumes  of  enemy  planes  or  missiles.  The  alternative  is  to  calculate 
the  appropriate  emissioiis  using  sophisticated  computer  codes.  A  number  of 
computer  codes  are  available  for  the  calculation  of  rocket  exhaust  plumes.  These 
computer  codes  have  not  been  fully  validated,  and  validation  of  the  codes  is 
required  to  increase  confidence  in  their  use  and  to  provide  data  for  their 
improvement.  Measurement  of  temperature  profiles  of  rcxrket  exhausts  offers  an 
important  method  for  validating  codes. 

Temperature  profiles  can  be  measured  using  intrusive  or  non-intrusive 
methods.  Non-intrusive  optical  methods  based  on  CARS  have  been  developed  in 
the  UK  for  temperature  measurement  of  particulate  free  phimes,  but  the 
methods  are  complex  and  expensive.  Intrusive  probes  based  on  thermocouples 
are  a  possible  alternative,  but  a  rocket  exhaust  presents  a  very  extreme 
environment  of  high  temperahire,  high  velocity,  high  turbulence  and  chemical 
reactivity.  These  factors  place  severe  constraints  on  the  type  and  design  of  probe 
that  can  be  used. 

In  the  current  work  an  intrusive  probe  was  the  only  prachcal  option.  There 
were  no  standard  thermocouple  probes  available  for  plume  measurement,  so  it 
was  necessary  to  develop  probes  in-house.  A  literature  survey  did  not  reveal  any 
suitable  designs,  but  much  ir\formation  was  foimd  on  which  to  base  a  design.  A 
twin  shielded  molybdenvun  probe  contaming  a  tugsten-rhenium  thermocouple 
was  designed  and  tested  with  a  small  scale  motor.  The  testing  of  the  probe  and  it 
use  in  measuring  temperature  in  the  plume  of  an  air-air  missile  are  discribed  in  a 
separate  paper(l]. 

This  paper  siunmarises  the  literahire  survey  on  high  temperature 
measurement  and  describes  the  design  of  probes  used  in  pliune  measurements. 
Corrections  for  conduction  losses  are  also  discussed. 


2.  Literature  Survey 


2.1  Thermocouple  recovery  and  correction  factors 

In  order  to  assess  the  results  of  experiments  on  thermocouple  behavioiu:,  it  is 
necessary  to  imderstand  the  concepts  of  recover^  factor  and  correction  factor.  A 
consistent  and  clear  discussion  is  given  by  Cambel  and  Jennings  in  [2],  and  a 
summary  of  their  discussion  is  given  below. 

The  deviation  of  the  measured  temperatiure  from  the  true  total  (or  stagnation) 
temperature  can  be  represented  by  a  correction  factor  K,  given  by 
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T 


Where  Tp  is  the  temperature  recorded  by  the  probe  which  may  exchange  heat 
with  the  envirorunent  through  radiation  and/or  conduction,  Tj  is  the  static 
stream  temperature  and  r^is  the  total  temperature.  If  the  calibration  factor  for  a 
probe  is  known,  as  well  as  the  heat  capacity  and  velocity  of  the  gas,  then  the 
total  temperature  can  be  calculated  from  a  measurement  of  the  probe 
temperatxue. 

If  it  is  assmned  that  the  temperatiure  ser\sor  exchanges  no  heat  with  its 
svuTOundings,  any  deviations  between  the  sensor  temperature  and  the  total 
temperature  of  the  flowing  material  will  be  due  to  the  flow  behaviour  in  the 
boundary  layer  around  the  ser\sor.  The  flow  impinging  on  the  front  of  the 
surface  of  the  sensor  is  brought  to  rest,  and  its  temperature  is  the  total,  or 
stagr\ation,  temperature,  T„,  which  is  given  by 

7;  =  7;  +  K-/2c, 


where  Cp  is  the  specific  heat  and  V  the  velocity.  However,  the  flow  aroimd  the 
remainder  of  the  sertsor  will  have  a  progressively  smaller  change  in  velocity 
compared  with  the  flow  at  the  tip,  and  the  local  temperatiue  will  be  closer  to  the 
stream  temperahire.  The  sensor  as  a  whole  will  therefore  recora  a  temperature 
which  is  between  Tj  and  because  the  flow  energy  is  only  partly 
recovered.  7^,^,  is  given  by 

7:^  =  r,+r.r^/2e, 

Reauranging  and  elinunating  /  2c^  gives 


The  factor  r  is  known  as  the  recovery  factor. 

If  a  thermometer  is  designed  such  that  heat  losses  are  negligible,  then  K  =  r  and 
=  Tp.  However,  this  is  an  approximation  and  K  is  not  synonymous  with  r. 
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Figure  1:  Dependence  of  various  temperatures  on  gas  velocit)/  from  Cambel  and 
Jenningsl2]. 


The  difference  between  the  two  factors  is  illustrated  in  figure  1  from  reference 
[2].  For  flow  with  constant  Ty ,  an  increase  in  velocity  decreases  the  static 
temperature,  and  the  adiabatic  sensor  temperature  Ty„,  decreases  by  a  lesser 
amount.  The  temperature  of  the  probe  with  heat  losses,  T.„  is  lower  than  T„y, .  At 
low  velocities  T^„  is  much  lower  than  T„j,,  because  heat  loss  dominates  the  heat 
gain  from  convection. 

As  noted  by  Moffat  [3],  the  performance  of  the  most  complex  probe  is 
determined  by  the  jimction  inside  it  and  the  environment  to  which  the  jimction 
is  exposed.  Radiation  shields  do  not  change  the  junction's  response  to  its 
environment,  they  change  the  environment. 

Moffat  [3]  reviewed  the  literature  on  the  recovery  factor  of  bare  wire 
thermocouples  and  recommended  values  of 

r  =  0.68  ±  0.07  wire  normal  to  flow 
r  =  0.86  ±  0.09  wire  parallel  to  flow. 

A  spherical  bead  much  larger  than  the  wire  diameter  moves  the  value  of  r 
towards  0.75  for  flow  in  both  directiorrs.  Rough  surfaced  wires  show  higher 
rccc'.’ery  factors  than  smooth  wires.  Moffat(3]  states  that  reproducibility  of 
recovery  factor  seems  to  be  about  ±3%  to  5%  for  junctions  made  with  normal 
care  by  experienced  technicians. 
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2.2  Performance  of  various  thermocouple  types 

Thermocouples  can  be  used  in  the  bare  wire  state  or  in  some  type  of  shield  to 
make  a  prc^.  An  imderstanding  of  the  behavior  of  bare  wire  thermocouples  is 
important  because  these  are  the  basis  of  shielded  probes. 


Bare  wire  thermocouples  were  studied  in  subsonic  flows  by  Scadron  and 
Warshawsky  [4].  The  thermocouples  used  are  illustrated  in  figiue  2  and  the 
dimensions  and  materials  are  given  in  Table  1.  Four  pairs  of  thermocouple 
nraterials  and  three  different  wire  diameters  were  used.  The  Mach  number  range 
was  0.1  to  0.9  and  the  Reynolds  munber  based  on  wire  diameter  ranged  between 
250  and  30,000.  The  gas  temperahire  was  573‘’K. 


Detailed  aruilyses  were  made  assiuning  that  bare  wire  thermocouples  were 
mounted  between  supports  with  temperatures  which  did  not  vary  significantly 
with  time.  Scadron  and  Warsawsky  made  a  large  number  of  measurements  of 
time  constants,  and  produced  detailed  calculations  for  evaluating  time  constants, 
conduction  corrections  and  radiation  correctiorrs.  However  the  subsonic 
corrections  they  derived  are  too  complex  for  routine  use. 

Configuration  I  Conligurotion  II 


Ceramic 

tube 


junction 


Figure  2:  Probes  from  Scadron  and  Warsltaiosky  [4]. 


Tabu  1; 


Probe 

Thermocouple  material 

Dimenilont  (mm) 
a  b  c 

d 

Configuration 

A 

Chromel-conslontan 

051 

.T.81 

1.27 

0.22 

I 

D 

•  • 

1.02 

7.62 

2..M 

0.45 

I 

C 

•  • 

254 

12.7 

655 

1.41 

U 

D 

Iron-oonatantan 

1.02 

7.62 

2.54 

0.41 

1 

E 

Chromel-alumel 

1.02 

7.62 

2.54 

0.39 

I 

F 

Pt-13%Ro/Pl 

1.02 

7.62 

2.54 

0.44 

I 

G 

Chromel<analantan 

• 

- 

- 

0.44 

Straight  wire 

A  summary  of  values  of  the  ratio  of  measured  temperature  to  total 
temperature,  and  static  tempera*,  ure  to  total  temperature,  versus  Mach  number 
is  given  in  figure  3. 


A  major  finding  was  that  over  die  range  of  conditions  studied,  the  Nusselt  | 

number  for  How  perpendicular  to  the  wire  was  given  by 


Nu  =  (0.427±0.018).Re<0‘'’‘'«>“!').PrOS 
or 

Nu  =  (0.478±0.002).R^O^Pr«■•'' 


where  the  variables  used  to  calculate  Re  are  evaluated  at  the  total  temjperature. 
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Figure  3;  The  ratio  of  measured  temperature  to  total  temperature,  arid  static 
temperature  to  total  temperature,  versus  Mach  number  from  Scadron  and  Viarshamda/ 
14J. 


Bare  wire  thermocouples,  as  well  as  shielded  probes,  were  analysed  by 
Stickney  [5].  Six  probes  were  tested,  and  they  are  illustrated  in  figure  4.  The  vent 
to  entrance  area  ratio  of  the  shielded  probe  6  was  0.43.  The  vent  in  probe  5 
consisted  of  2  unused  holes  in  the  insulator  for  thermocouple  lead  wires  which 
were  not  dimensioK&d.  The  wire  diameter  was  0.305  mm,  and  assuming  that  the 
hole  waf  0  4  mm  in  diameter,  the  area  ratio  would  be  0.2.  Probe  4  was  only 
partially  shielded,  so  a  vent  to  entrance  ratio  is  not  relevant. 

So  called  recovery  factors  (r),  which  are  referred  to  as  calibration  factors  (K) 
here,  were  measured  over  the  Mach  number  range  0.2  to  2.2  at  a  total 
temperature  of  21-38°C.  At  these  temperatures  heat  flow  from  the  probe  would 
be  negligible,  and  r  can  be  equated  to  K  and  measured  directly.  The  measured 
recovery  factors  are  plotted  in  figure  5. 


Figure  4(a);  Probes  1  and  2  from  Stickiie]/  15J. 


Figure  4(b):  Probe  3  from  SHckitey  15]. 
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Figure  4(e):  Probe  6  from  Stickney  [SJ. 


From  figure  5  it  can  be  seen  that  the  unshielded  probes  were  strongly  affected 
by  the  transition  from  subsonic  to  supersonic  flow.  Probe  2  shows  that  while  the  ^ 

probe  remains  behind  the  bow  shock  the  curve  is  a  continuation  of  the  subsonic 
curve,  but  when  the  sensor  breaks  through  there  is  a  sudden  drop  in  recovery 
ratio.  The  sensors  in  the  shielded  probes  always  remain  behind  the  bow  shock. 

The  shielded  probes  showed  smaller  variations  of  recovery  factor  with  probe 
design,  and  the  double  shielded  probe  5  showed  the  least  effect  of  changing 
Mach  number.  * 

These  results  show  that  for  flows  involving  both  supersonic  and  subsonic 
regions  shielded  probes  must  be  used. 

A  low  temperature  shielded  probe  for  supersonic  flows  was  developed  by  • 

Goldstein  and  Scherrer  [6].  A  section  through  the  probe  is  illustrated  in  figure  6. 

The  main  part  of  the  probe  was  made  from  lucite,  and  the  sensor  was  placed  just 
above  the  hemispherical  surface  at  the  back  of  the  probe  because  this  was  where 
the  flow  was  nearly  stagnant,  but  there  was  sufficient  flow  for  heat  transfer  to 

the  sensor.  ^ 

The  air  sample  enters  the  blunt  shield  at  a  subsonic  velocity  behind  the 
detached  bow  wave.  It  is  necessary  to  have  a  high  Reynolds  number,  based  on 
entrance  diameter,  so  that  the  boimdary  layer  does  not  fill  the  interior  of  the 
probe.  If  the  boundary  layer  approaches  the  duct  radius  the  conduction  to  the 

shield  will  be  felt  by  the  serrsor.  * 


Mach  number 


Figure  5:  Therwocouple  calibration  factors  from  [5]. 


Goldstein  and  Scherrer  listed  two  problems  in  a  shielded  probe:  1)  getting  a 
fluid  sample  to  the  thermocouple  with  the  minimiun  heat  loss  through  the 
shield,  and  2)  arranging  the  thermocouple  sensor  in  the  probe  so  that  it  will 
receive  the  maximum  amount  of  heat  from  the  sample  and  loose  the  minimum 
amoimt  of  heat  through  radiation  and  conduction. 

The  total  temperature  during  the  tests  was  approximately  15.6‘^'C  and  the  Mach 
number  was  varied  between  1.2  and  2.1.  The  calibration  factor  was  found  to  be 
0.992  over  the  Mach  number  range  1.36  to  2.01.  The  optimum  ratio  of  vent  area 
to  entrance  area  was  between  0.5  and  0.625  for  a  Mach  number  of  15. 

A  range  of  probes  were  tested  in  subsonic  flows  by  Glawe,  Simmons  and 
Stickney  (7].  All  sensors  were  chromeTalumel  thermocouples.  Nine  types  of 
probe  were  tested.  Probe  1  is  illustrated  in  figure  7(a).  Probe  2  will  not  be 
considered  here.  Probes  3  to  6  were  similar  to  probes  tested  by  Stickney  [5]:- 

Probe  3  similar  to  Stickney  probe  4,  figure  4(c) 

Probe  4  similar  to  Stickney  probe  1,  figure  4(a) 

Probe  5  similar  to  Stickney  probe  3,  figure  4(b) 

Probe  6  similar  to  Stickney  probe  6,  figure  4(e) 
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Figure  7(a);  Probe  1  from  Glawe,  Simmons  and  Slickney  [7], 
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Figure  7(b):  Prolv  7  from  Glawe,  Simmom  and  Stickney  [7], 
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Figure  7(c):  Probe  8  from  Glawe,  Simmons  and  Stickney  17], 
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Probes  7  and  8  are  illustrated  in  figures  7(b)  and  7(c).  Probe  9  was  similar  to 
probe  8,  but  it  was  slightly  larger.  Probes  1,  8,  and  9  were  aspirated,  and  it  was 
found  that  a  constant  indicated  temperatiue  of  probe  9  was  obtained  when  the 
pressure  across  the  probe  was  greater  than  100  mm  of  Hg. 

Temperatiue  and  time  constant  data  were  obtained  at  static  pressures  2/3  to 
4/3  atmospheres,  temperatures  833  to  1389°K,  and  Mach  numbers  0.3  to  0.9. 
Recovery  factors,  defined  by  A  =  (Tt  -  Tg)/Tt,  were  obtained  for  Mach  numbers 
from  0.2  to  0.9  and  static  pressures  from  1/7  to  4/3  atmospheres  and  at  room 
temperatiue. 

True  total  temperatures  were  assumed  by  Glawe,  Simmons  and  Stickney  [7]  to 
be  giv'en  by  probe  1,  and  the  deviation  of  the  temperatures  given  by  the  other 
probes  from  the  true  total  temperature  were  derived  from  comparisons  with 
probe  1.  Temperatiue  deviations  at  a  Mach  number  of  0.3  are  given  in  figure  8. 
Probes  1,8  and  9  showed  no  variation  of  recovery  factor  with  Mach  number  over 
the  range  0.2  to  0.9,  but  all  the  other  probes  showed  considerable  variation. 

The  transient  response  for  probes  8  and  9  strongly  deviated  from  a  first  order 
system.  The  temperatiue  initially  rose  at  a  high  rate  but  the  rate  dropped  as  the 
final  temperatiue  was  approached.  It  was  suggested  that  the  later  part  of  the 
response  was  due  to  the  slower  heating  of  the  shields  w’hich  are  of  considerably 
greater  mass  than  the  thermocouple. 

Glawe,  Simmons  and  Stickney(7]  derived  an  expression  for  the  radiation 
correction  (including  a  conduction  correction)  and  for  the  time  constant. 

Design  details  and  test  results  were  given  for  two  types  of  total  temperature 
probes  for  hypersonic  boimdary  layer  measurements  were  given  by  Albertson 
and  Bauserman(8].  The  boundary  layer  to  be  studied  was  approximately  25  mm 
thick,  so  a  small  probe  was  needed.  Their  probe  B  is  illustrated  in  figure  9.  Probe 
A  was  similar  to  B  with  a  support  sleeve  around  the  outside.  The  outside 
diameter  of  A  was  2.4  mm  and  B  was  1.5  mm.  To  reduce  vertical  averaging  the 
entrance  to  the  probes  was  flattened,  probe  A  was  0.75  mm  high  and  probe  B 
was  0.5  mm  high.  The  thermocouple  was  platinum-13%-  rhodium  vs  platinum 
(type  R),  wire  diameter  0.25  mm.  A  single  platiniun-20';^.-rhodium  shield  0.125 
mm  thick  was  used  in  both  cases,  the  main  difference  was  the  vent  to  entrance 
area:  A  50%  and  B  60%.  Estimated  Reynolds  number  in  the  probe  ranged  from 
63  to  2470. 
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Figure  8:  Temperature  corrections  at  Mach  number  03  for  probes  Glaroe,  Simmons 
and  Stichiey  [7] 
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Figure  9:  Albertson  and  Bauserman  probe  B 18J. 


The  tests  were  carried  out  at  the  boundary  layer  edge  where  the  Mach 
niunbers  were  5  and  6.2  and  near  die  wall  where  the  Mach  number  was  1.7.  The 
pressures  were  equivalent  to  altitudes  ftom  80,000  to  126,000  ft  and  the  nominal 
stream  temperatures  were  1444  K  and  1777  K. 


Albertson  and  BausermanlS)  used  a  suggestion  by  Wmkler(9]  that  the 
calibration  foctor  K  is  proper titxial  to  Nusselt  munber  given  by:- 

where  lip  is  the  convective  heat  transfe'  coefhciait  to  the  junction,  dy,  is  the 
junctian  wire  diameter,  ai\d  Ic^  is  the  conductivity  of  the  wire.  Winkler  also 
suggested  a  correlation  of  K  with  PtTj'^  where  p^  is  the  total  pressure  in  die 
probe  (They  plotted  die  relation  but  there  was  a  large  amount  of  scatter,  and 
die  relation  was  not  convincingly  established  in  this  case.) 

Albtttson  and  Bauserman(8]  derived  a  conduction  correction,  and  more  details 
of  it  will  be  given  in  Section  4.2  because  it  is  the  basis  for  the  corrections 
developed  in  the  present  work.  The  conduction  correctitm  is  a  fimction  of 
lengdi/ diameter  of  the  thermocouple  wires,  the  thermophysical  properties  of  die 
wires,  convective  heat  transfer  to  ^  wires  and  the  probe  support  temperahure. 

They  also  calculated  a  radiation  correction,  which  is  a  function  of  the 
emmitance  of  the  thermocouple,  the  shield,  the  radiation  view  factor  between  the 
shield  and  the  wires,  and  the  temperature  of  the  radiation  shield. 

A  water  cooled  probe  for  supersonic  flows  was  developed  by  Lagen  and  Seiner 
(101  and  it  is  illustrated  in  figure  10.  The  probe  was  made  from  AISI  Type  347 
stainless  steel,  and  it  was  water  cooled  and  mounted  on  a  water  cooled  wing  to 
increase  the  temperatures  the  probe  coidd  sustain.  The  mass  flow  rate  of  coolant 
in  die  probe  was  0.0295  kg/s.  The  thermocouples  were  type  K  chromel-alumel. 

Total  temperature  measurements  were  made  tor  liow  from  a  89.9  mm  diameter 
nozzle  with  a  Mach  number  of  2  and  temperatiures  up  to  1366K. 

The  surface  temperature  of  the  probe  in  operation  was  measured  with  a  high 
resolution  infrared  system.  The  emissivity  of  the  probe  surface  was  foimd  to  be 

om 

Probes  were  placed  in  the  flow  for  tests  of  duration  of  10  seconds  with  both  the 
probe  and  wing  with  coolant  dxculating  and  not.  Plots  of  measured  temperatiue 
vs  total  temperature  are  given  in  figure  11.  Radiaticm  and  conduction  factors 
calculated,  and  die  corrected  results  are  also  plotted  in  figure  11.  The  cooled 
probes  reached  steady  state  in  7-10  seconds,  but  separate  tests  showed  that  the 
uncooled  probes  took  up  to  2  minutes  to  reach  steady  state. 


14 


Figure  10:  Probe  from  Lagen  and  Seiner  [lOJ. 


Figure  11:  Measured  or  calculated  temperature  versus  true  total  temperature  from 
Lagen  and  SeinerllO}.  0  -  Calibration  points,  A  -  measured  bead  temperature,  []  - 
calculated  temperature, - ideal  calculated  temperature. 
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Figure  12:  Bartlett,  Edivards  and  Hillier  [11]  probes. 


A  total  temperature  probe  for  a  gun  tiumel  was  developed  by  Bartlett, 
Edwards  and  Hillier  (11].  A  schematic  of  the  probe  is  given  in  figure  12.  The 
bask  design  is  a  fine  thermocouple  stretched  between  2  support  prongs.  Total 
temperature  measurements  were  made  in  a  hypersonic  boundary  layer  with  run 
times  of  order  3-4  ms.  They  quote  a  formula  for  correction  including 
compressibility  effects  assuming  no  loss  to  surroimdings.  The  calculated 
corrections  were  considered  too  idealised  and  experimental  calibration  was 
required. 

The  free  stream  conditions  were  To  =  1064®K,  =  9.26,  Re^  =  5.5x10^. 

The  first  probe  is  illustrated  in  figure  (12a),  but  the  recovery  factor  was  too  low. 
It  was  improved  by  iiKreasmg  the  L/D  and  moving  the  bead  away  from  the 
surface  of  the  insulated  support,  see  figure  (12b).  This  produced  a  great 
iiiq>rovement,  but  the  best  design  was  found  to  be  figure  (12c),  where  fire  L/D 
ratio  of  tt>e  thermocouple  was  180/1.  The  entrance  was  flattened  to  improve 
resolution  in  the  vertical  directiotv 

Moeller,  Noland  and  Rhodes  [12]  describe  a  number  of  probes  devdoped  at  or 
for  NASA  and  NACA.  A  design  for  a  foermocouple  to  operate  at  >3250  K,  see 
figure  13,  resulted  horn  a  contract.  It  consisted  of  a  W  sl^th  and  a  coaxial  W- 
26Re  centre  ccmductor.  The  sheath  was  formed  by  a  thermochemical  process.  A 


mandrel  was  made  from  mild  steel  and  the  W-26Re  wire  was  placed  in  the 
centre  longihidinal  hole  and  crimped  at  the  end.  The  assembly  was  placed  in  a 
quartz  tube  and  heated  by  induction.  Vaporised  timgsten  hexafluoride  was 
introduced  as  the  mandrel  was  rotated.  When  the  coating  was  thick  enough  the 
assembly  was  placed  in  hot  concentrated  hydrochloric  acid  to  dissolve  the 
mandrel.  However  use  of  this  probe  was  restricted  because  it  was  susceptible  to 
stress  fatigue. 
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Figure  13:  Probe  from  Moeller.  Noland  and  Rhodes  112].  % 

No  high  temperature  sealing  materials  were  available,  so  the  thermocouple 
wires  were  brought  out  to  a  low  temperature  region  where  conventional 

materials  could  be  used.  Compensated  lead  wires  were  tested,  but  they  ^ 

introduced  errors  unless  they  were  used  at  temperatures  less  than  about  370°K. 

Preheating  the  thermocoviples  to  near  the  expected  test  temperature  was  foimd 
to  decrease  the  response  time. 

Two  types  of  shielded  thermocouples,  see  figure  14,  were  used  in  a  wind  • 

hmnel  with  total  temperature  of  1820"K  and  a  static  pressure  range  from  15  to 
120  microns  of  Hg.  The  thermocouples  and  shields  were  made  from  Pt-lORh/Pt. 

The  cold  shield  probe  minimised  conduction  losses  by  using  small  diameter 
wires  supported  on  V  shaped  thicker  thermocouple  wires.  The  shield  had 

sufficient  mass  not  to  heat  up  during  the  run,  so  radiation  corrections  could  be  0 

calculated  with  fixed  boimdary  conditions.  A  probe  with  3  radiation  shields  is 
shown  in  figure  14. 
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A  fast  response  probe  for  measuring  temperatures  up  to  1760^K  for  10  seconds 
in  gases  with  velocihes  up  to  240  m/s  and  pressures  of  5.6  MPa  is  illustrated  in 
figiure  15.  The  hme  constant  was  0.1  second.  A  low  thermal  conductivity  alloy 
was  used  for  the  moimting  stem  to  reduce  heat  loss  to  the  walls.  Pt-lORh 
thermocouples  were  used  because  W/W-Re  thermocouples  were  too  brittle  and 
the  hmgsten  was  too  hard  to  fabricate. 
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Figure  14:  Tivo  radiation  shielded  thermocouple  probes  from  Moeller,  Noland  and 
Rhodes  112]. 
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Figure  15:  A  fuat  n'^poiiJiiig  >liichicJ  i>robc/roiu  Moeller,  NolauJ,  aiul  RhodefllS], 


2.3  Probe  materials 

Moeller,  Noland  and  Rhodes!  12]  reviewed  the  experience  in  NASA  in 
temperature  measurement  in  areas  such  as  cryogenic  temperatures,  high 
temperatures,  gas  temperahires,  surface  temperature  and  instrumentation 
techniques.  Some  of  the  points  of  interest  here  are  summarised  below. 

The  noble  metal  alloy  platinum-rhodium(Pt-Rh)  thermocouples  canjiot  be  used 
much  above  192(h’K,  and  iridiiun-rheniiun/iridium(Ir-Rh/Ir)  thermocouples  are 
limited  to  temperahires  below  242(h’K.  Timgsten-Rhenium  thermocouples,  either 
W/W-26Re  or  W-5/W-26Re,  can  be  used  up  to  3255'^’K  but  the  W-5/W-26Re  is 
preferred  because  it  is  less  brittle.  These  thermocouples  have  high  emf,  chemical 
stability  at  high  temperatures  in  vacuum  and  inert  atmospheres,  and  relatively 
low  vapour  pressures. 

Tables  of  physical  and  thermoelectric  properties  of  high  temperature 
thermocouple  materials  are  given  by  Caldwell  |13].  The  use  of  refractory  metals 
for  ultra  high  temperature  thermixouples  is  discussed  by  Lachman  and 
McGurty  (14].  In  particular,  data  on  W/Re  and  W/Mo  thermocouples  are 
given.  W/Mo  thermocouples  have  the  disadvantage  that  the  emf  vs 


temperature  curves  have  a  minimum  value  near  1000  K  and  the  thermal  emf  is 
low. 

There  are  a  number  of  high  temperature  insulation  materials  that  can  be  used 
in  probes.  Moeller,  Noland,  and  Rhodes  [12]  consider  beryllia  (BeO)  melting  at 
2820"K,  and  magnesia  (MgO)  melting  at  3070‘’K.  The  electrical  resistance  falb 
dramatically  with  increasing  temperature.  In  a  particular  geometry  the  resistance 
in  a  magnesia  probe  fell  from  lO'^  ohms  at  lOOO^K  to  200  ohms  at  2200°K. 
Molybdeniun  reacted  with  BeO  at  2600°K,  but  tungsten  did  not.  At  about  2700°K 
timgsten,  tantalum,  and  W-26Re  wires  react  slightly  with  BeO. 


3.  Thermocouple  Probe  Design  Considerations 

Moffat  [3]  stated  that "  the  design  of  a  probe  for  a  given  acairacy  is,  in  reality, 
the  design  of  a  device  to  provide  an  environment  inside  the  probe  which  will 
allow  a  bare  wire  junction  to  measure  gas  temperature  with  acceptable 
accuracy". 

The  two  main  sources  of  error  in  the  measurement  of  total  temperature  are 
conduction  and  radiation  from  the  thermocouple.  As  well  as  conduction  of  heal 
from  the  junction  to  the  probe  support  base,  Moffat  stated  that  the  slueld  and 
support  paiallel  to  the  flow  have  a  recovery  factor  of  0.86,  the  same  as  a  bare 
junction  exposed  to  the  stream.  This  depression  of  temperature  below  the  gas 
total  temperature  provides  an  additional  driving  force  for  conduction  and 
radiation  loss  from  the  junction. 

Tne  conduction  error  can  be  reduced  by  (3]: 

1.  Decreasing  the  difference  in  temperature  betw-een  the  junction  and  the 
moimt. 

2.  Increasing  the  junction  length. 

3.  Decreasing  the  wire  diameter  or  thermal  conductivity. 

4.  Increasing  the  heat  transfer  coefficient  h^. 

Only  the  junction  length  can  easily  be  changed.  The  heat  transfer  coefficient  can 
be  increased  by  increasing  the  internal  gas  velocity  in  the  probe,  but  this  can 
only  be  carried  so  far  before  the  flow  can  no  longer  be  considered  stagnant. 

The  radiation  error  can  be  reduced  by  [3]: 

1.  Decreasing  emissivity. 

2.  Increasing  h^. 

3.  Increasing  the  wall  temperature. 

The  emissivity  can  be  decreased  by  coating  the  junction  with  a  material  hke 
platinum,  but  contamination  can  increase  the  emissivity.  h^  can  be  increased 
above  the  value  obtained  from  the  free  stream  pressure  alone  ^v  .aspirating  the 


probe  using  an  external  piunp(3].  The  wall  temperature  can  be  increased  by 
using  multiple  shields. 


From  the  literature  reviewed  above,  it  is  apparent  that  bare  wire 
thermocouples  are  suitable  or\ly  for  low  temperature,  subsonic  flows.  Single 
shields  have  been  used,  but  at  high  temperatures  the  radiation  losses  wouid  be 
too  high  or  calculated  correction  factors  imcertain.  The  multiple  shield  type  of 
figure  14  would  be  an  improvement,  but  the  shields  would  only  equilibrate  to 
the  recovery  temperatiure  and  not  the  total  temperature,  and  at  high 
temperatures  this  difference  could  be  significant.  The  aspirated  probes  appeared 
to  give  acceptable  performance,  but  the  difficulty  of  providing  aspiration  rules 
them  out  for  the  present  purpose.  The  side  entrance  multishield  probe  in  figure 
15  appeared  to  be  the  best  of  the  probes  surveyed,  but  it  would  be  too  difficult  to 
manufactiure  from  timgsten,  which  would  be  required  at  the  temperatures  in 
rocket  exhausts. 
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Figure  16:  Schematic  of  probe  deveLped  for  tin  ■■  toork. 


The  final  design  chosen  combined  the  requirements  of  multiple  shields, 
sturdiness  and  practicality  of  manufacture.  The  probe  is  illustrated  in  figure  16. 
It  consists  of  a  twisted  wire  W-5%Re/w-26%Re  thermocouple  supported  by  a 
twin  bore  alumina  rod.  The  thermocouple  wire  diameter  is  0.5  mm  and  the 
shield  material  was  tungsten. 

For  the  purposes  of  calculating  the  conduction  from  the  probe  it  would  have 
been  preferable  to  have  used  longer,  straight  rather  than  twisted,  thermocouple 
wires  in  the  probe.  However,  there  was  concern  about  the  strength  of  that 
configuration  and  it  was  decided  to  use  the  twisted  wires.  Subsequent 
experience  has  indicated  that  straight  wires  may  have  been  adequate. 

The  irmer  shield  is  vented  with  4  holes  with  a  total  area  of  one  half  the  inner 
shield  entrance  area.  The  space  between  the  shields  forms  a  second  chamber 
which  is  vented  through  4  holes  with  an  area  of  half  the  entrance.  The  gas  in  the 
second  chamber  is  slowed  and  heated  to  near  the  temperatiue  of  the  inner 
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chamber.  Observatioi\  of  the  irmer  shield  through  an  outer  vent  hole  during  a 
test  confirmed  that  the  inner  shield  heated  to  incandescence  very  quickly, 
confirming  that  the  shield  was  at  a  temperature  high  enough  to  minimise 
radiation  losses  from  the  thermocouple  jimction. 

There  was  concern  that  the  alumina  support  rod  would  not  provide  sufficient 
electrical  imulation  imder  measurement  conditions.  Two  unconnected 
thermocouple  wires  were  inserted  in  a  rod  and  the  rod  heated  imder  an  oxy- 
acetylene  flame  for  2  minutes.  The  electrical  resistance  between  the  wires  did  not 
drop  below  8  megohms,  so  the  performance  of  the  alumina  was  taken  to  be 
satisfactory. 


4.  Conduction  and  Radiation  Corrections 

4.1  General 

The  temperature  indicated  by  a  probe  results  from  a  heat  flow  balance 


Qc  +  Qr  +  Qk  +  Qs  =  0, 

where 

Qj.  is  the  rate  of  heat  convection  between  gas  and  sensor 
Q,  is  the  rate  of  radiant  heat  exchange  between  the  sensor  and 
surroundings. 

Qk  is  the  rate  of  conduction  of  heat  from  the  sensor  by  the  thermocouple 
wire. 

Qg  is  the  rate  of  heat  storage  in  the  sensor. 

The  heat  convection  is  given  by 

Qc  -  he  (Tg  -  T^v), 

where  1\.  is  the  convective  heat  transfer  coefficient,  Tg  and  Ty^  are  the  gas  and 
junction  temperatures. 

Moffat  [3]  stated  that  convection  of  heat  to  the  thermocouple  from  the  gas 
stream  flowing  over  a  solid  surface  is  affected  by  the  boundary  layer  which 
becomes  established  and  it  prevents  the  free  stream  from  coming  in  contact  with 
the  surface.  The  thermal  conductivity  of  the  boundary  layer  determines  the  rate 
of  transfer  of  heat  from  the  free  stream  to  the  surface.  This  conductance  is 
described  by  the  heat  transfer  coefficient  1\. ,  which  is  a  fimction  of  Reynolds 
number  Re,  fluid  properties,  geometry  and  stream  turbulence.  The  Nusselt 
number  is  defined  as  Nu  =  h^  d  /  kf,  where  kf  is  the  thermal  conductivity  of 
the  fluid.  The  Prandtl  number  is  defined  as  CpP/kf ,  and  was  taken  to  be  equal  to 
0.7.  Moffat  13]  presented  the  following  empirical  relationships  between  Nu  and 
Re  for  Re  in  ^e  range  100  to  10000. 


Wire  normal  to  flow:  Nu  =  (0.44  ±0.06)  Re*’ “* 

Wire  parallel  to  flow:  Nu  =  (0.085  ±  0.009)  Re*’ 

The  data  for  these  empirical  relationships  were  obtained  at  relatively  low 
turbulence  levels,  and  rocket  motors  would  have  higher  turbulence  and  hence 
higher  Nusselt  numbers.  However,  use  of  the  empirical  relationships  would  lead 
to  a  conservative  probe  design,  and  this  was  considered  to  be  acceptable  in  the 
present  case. 

The  radiative  heat  trar\sfer  is  given  by 

where  a  is  the  Stefan-Boltzmann  corrstant,  is  the  emissivity  of  the  wire  and 
Tj  is  the  wall  temperature.  With  the  probe  design  considered  here,  it  was 
considered  that  a  radiation  correction  would  be  too  difficult  to  develop,  and 
since  the  probe  was  designed  to  minimise  radiation  losses,  a  radiahon  correction 
was  not  corwidered.  The  neglect  of  a  radiation  correction  would  increase  the 
uncertainty  of  the  results  by  an  unknown,  but  probably  not  significant,  amount. 

Heat  storage  is  given  by 
Qs  ~  Pw*-W  ^Ty^/9t 

where  is  the  density  of  the  wire  and  c^^  is  the  specific  heat  of  the  wire.  Under 
steady  state  conditions  heat  storage  effects  would  be  zero,  and  non-steady  effects 
will  not  be  considered  here. 

Heat  conduction  will  be  discussed  fully  in  the  next  section. 

4.2  Calculation  of  the  Conduction  Correction 

The  calculation  of  the  conduction  correction  is  based  on  the  method  described  by 
Albertson  and  Bauserman  [8].  Flow  within  the  probe  is  considered  to  be 
isentropic.  The  entry  region  flow  is  modelled  as  a  divergent  nozzle,  with  no  heat 
loss  in  the  recircxilating  flow  behind  the  expansion.  The  aft  end  is  modelled  as  a 
convergent  nozzle  with  an  exit  area  equal  the  sum  of  the  areas  of  the  vent  holes. 
The  flow  is  calculated  assuming  that  the  static  pressure  in  the  vent  is  equal  to  the 
external  stream  pressure.  An  equivalent  assmnption  that  the  velocity  through 
the  vent  holes  may  be  taken  as  equal  to  the  free  stream  velocity  if  the  axis  of  the 
tube  is  parallel  to  the  main  stream  and  there  are  no  obstructions  outside  the 
holes  was  suggested  by  Moffat  [3]. 

The  thermocouple  in  the  probe  will  be  modelled  by  the  idealised  geometry 
illustrated  in  figure  17. 
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Insulation 


Figure  17:  Schematic  of  probe  for  the  coriduction  correction. 


The  quantities  required  to  evaluate  the  correction  are:- 

hp  heat  transfer  coefficient  to  the  thermocouple 
Cp  specific  heat  of  the  gases  flowing  over  the  probe 
k„  thermal  conductivity  of  the  thermocouple  wire 
k-  thermal  conductivity  of  the  gas  in  the  probe 
A  area  of  the  thermocouple  wire  =  1.964x10'^  m^ 
diameter  of  the  thermocouple  wire  =  0.5x10*^  m. 
p  gas  density 

U  gas  velocity 

L  length  of  thermocouple  wire  diameter 
gas  viscosity 

y  ratio  of  specific  heats 

R  gas  constant 

relevant  temperatures 

and  the  dimensionless  numbers 

I 

Bi  Biot  number  /  4^,  )- 

Nusselt  number  h^d^  / 

Rf  Reynolds  number  pl/L  /  p 
Py  Prandtl  number  pcp/k^ 

The  heat  flow  in  a  small  segment  of  wire  was  assumed  to  be  given  by  a  one 
dimensional  heat  balance:- 


V,[7’a.-W]=-M[0] 


where  is  the  temperatiure  of  the  wire  if  there  is  no  heat  loss  to  the 

enviroiunent.  A  general  solution  of  the  equation  is 
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where  Bi  is  the  Biol  number  given  by  ( /  4A.  )- .  If  the  temperature  at  n  =  Lj 
was  known  the  conduction  correction  f -  1)  could  be  evaluated.  Since  this 
temperature  was  not  measured,  heat  traiisfer  from  n  =  L]  to  n  =  L2  was 
calculated  by  assuming  one  dimensional  heat  transfer  with  a  temperatiure  at  L2 
of  27^.  With  these  assumptions  the  correction  becomes 


{T^-L)  = 


r^-T{n  =  L,) 
A 


where  -i-  D^{L.  -  L^) 


A  =  cosh 


i  \-  {  4/?„  V 


(JJ. 


sinh 


4^ 


+  B.  cosh 
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(4) 

(5) 
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For  gas  flow  parallel  to  the  thermocouple  w'ires,  Moffat  (3)  demonstrated  a 
correlation  between  Nusselt  number  and  Reynolds  number  as 

A  =0. 095/^;' (7) 

The  quantities  h^, ,  c^, ,  lc„, ,  p  ,  and  p  depend  on  the  gas  composition  and 
temperature.  The  gas  composition,  c^„  gas  molecular  weight,  tind  velocity  U  are 
obtained  from  the  output  of  the  REP  code.  Ideal  gas  behaviour  is  assumed,  and 
the  value  of  y  is  obtained  from  c^,  and  R,  the  gas  constant.  The  gas  density  p  is 
calculated  from  the  pressure  using  the  ideal  gas  equation.  The  static  pressure  in 
the  stream  is  assumed  to  be  atmospheric.  The  convective  heat  transfer  coefficient 
hp  is  calculated  from  the  Nusselt  number. 

The  major  gases  in  the  plume  are  N2,  O2,  CO2,  H2O,  CO  and  HCl. 
Thermodynamic  data  for  HCl  was  not  readily  available,  and  it  was  a  relatively 
small  component,  so  HCl  was  not  included  in  calailations  of  gas  properties. 
Data  for  CO  was  available  only  over  a  very  limited  range,  and  since  it  was  a 
relatively  small  component,  it  too  was  not  included.  In  view  of  the  limited  data 
available  over  the  temperature  range  of  interest,  it  was  decided  to  approximate 
the  gas  as  a  mixture  of  air,  H2O  and  CO2  with  the  appropriate  composition. 
Literature  values  for  viscosity  of  air  were  obtained  for  temperahires  up  to 
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1873°K,  H2O  up  to  773'^K,  and  CO2  up  to  1325“K  (15].  Literatxire  values  of 
thermal  conductivity  data  for  air  and  H2O  were  obtained  for  the  same 
temperatures,  and  for  CO2  to  819‘’K  [15]. 

The  viscosity  and  thermal  conductivity  data  were  fitted  by  a  linear  hmcbon 
over  the  temperatiure  range  of  interest,  extrapolating  where  necessary.  The 
resulting  functions  are  given  below. 


Gas 

Viscosity 

Thermal  conductivity 

(I'as) 

(W/mK) 

Air 

2.0xUr‘'  +  2.25xU)-><T 

0.02  +  4..1.3xUr^ 

H2O 

.3.6xl(r«T 

-0.04  +  1.4.SxlO“*T 

CO2 

1..SX10-'' +  2.4xU)-^T 

0.(K)5  +  6.3.SxU)-^  T 

The  thermal  conductivity  of  the  thermocouple  was  taken  to  be  that  of  tungsten. 
The  conductivity  values  obtained  from  Raznjevic{15]  show  a  decrease  with 
increasing  temperatiue  to  a  value  of  98.9  Wm'^K'^  at  1273"K,  and  then  a  nearly 
linear  increase  to  146.5  Wm'^K'l  at  2673‘'K.  However,  Deshpande  and 
Taylor(16]  report  the  opposite  trend  at  high  temperatures,  with  the  conductivity 
dropping  from  116  W/mK  at  1200‘’K  to  about  80Wm"^K'^  at  3000'^’K.  The  values 
of  Raznjevic  were  chosen  because  these  would  give  greater  calculated  corrections 
cind  reflect  the  greater  imcertainty  in  the  temperatiue  value.  The  high 
temperature  branch  was  approximated  by 

=  55  +  3.5x10-2t  Wm'^K-l 

For  a  particular  temperature  in  the  probe  the  material  parameters  were 
evaluated  and  the  Reynolds  and  Prandtl  numbers  were  calculated.  The  Nusselt 
number  was  calculated  from  equation  7,  and  this  value  was  used  to  calculate  the 
convective  heat  transfer  coefficient  /i,,.  The  Biot  number  and  the  factors  D3,  D4 
and  D5  were  then  calculated  and  used  to  evaluate  the  conduction  correction 
from  equation  3. 

The  method  used  to  calculate  the  temperature  and  pressure  in  the  probe 
depended  on  whether  the  stream  was  supersonic  or  subsonic.  These  two  cases 
wiU  be  considered  below. 


Supersonic  stream 

When  the  external  stream  flow  is  supersonic  there  is  a  detached  shock  before  the 
probe  which  can  be  considered  to  be  a  normal  shock  in  front  of  the  probe 
entrance.  The  flow  variables  at  various  locations  will  be  identified  by  the 
following  subscripts 

1  Upstream  of  the  shock 

2  Downstream  of  the  shock,  outside  the  probe 

3  Inside  the  probe 


where  the  subscript  v  refers  to  the  conditions  in  the  vent.  The  vent  is  relatively 
close  to  the  entrance  to  the  probe,  so  the  flow  past  the  vent  will  be  similar  to  the 


flow  behind  the  shock,  and  hence  M,,  =  M2.  The  area  ratio  in  the  current  case 
equals  2,  so  the  Mach  number  in  the  probe,  M ),  can  be  evaluated. 

The  velocity  in  the  probe  is  given  by 

U,=M,ylyRT;  (14) 

Once  U3  has  been  determined,  the  conduction  correction  can  be  calculated 
from  equation  3.  Since  the  flow  in  the  probe  is  not  stagnant,  the  corrected 
junction  temperature  is  not  the  total  temperature  but  the  recovery  temperatiue. 
A  factor  (1  -  r)Uy~  f 2c ^  must  be  added  to  the  measured  temperature. 

The  corrections  were  calculated  using  a  spreadsheet,  and  a  typical  output  is 
illustrated  in  figure  18.  Pj,  Uj,  gas  composition,  and  other  gas  properties  were 
obtained  from  REP  output.  Tj  had  to  be  calculated  iteratively  since  the  value  of 
Mj  was  required  in  equation  8,  and  a  value  of  Tj  was  needed  to  calculate  Mj 
from  equation  10.  Trial  values  of  Tj  were  substituted  into  equations  8  and  9  until 
the  correct  value  of  Tyj  was  obtained. 

The  values  of  the  variables  after  the  shock  in  region  2  were  calculated  directly. 
The  value  of  M3  in  region  3  of  the  probe  w-as  obtained  by  trial  and  error  by 
calculating  the  entrance/vent  area  ratio,  and  adjusting  M3  until  the  area  ratio 
was  the  correct  value.  The  calculation  of  the  corrections  was  then 
straightforward. 

Subsoitic  stream 

In  this  case  region  1  does  not  occur,  and  the  variables  in  region  2  are  calculated 
in  the  same  way  as  the  variables  in  region  1  in  the  supersonic  stream  case.  The 
variables  in  region  3  in  the  probe  are  c^culated  in  the  same  way  as  the  variables 
in  region  3  in  the  supersonic  stream. 

4.3  Accuracy  of  probe  temperature  measurements 

It  would  be  desiiable  to  be  able  to  estimate  the  accuracy  of  the  temperatures 
measured  by  probes  described  here.  However,  it  is  not  practical  to  provide  such 
an  estimate.  TTie  sources  of  deviation  of  the  measured  value  from  the  true  value 
have  been  analysed  and  shown  to  be  complex.  For  the  conduction  correction  the 
values  of  the  physical  properties  were  in  some  doubt,  and  there  would  be 
variability  in  the  construction  of  the  probe.  The  effect  of  these  uncertainties  on 
the  magiutude  of  the  correction  cannot  be  calculated.  Losses  due  to  radiation  are 
not  known,  but  are  mirumised  by  the  probe  design. 

The  assumptions  made  in  the  development  of  the  probe  design  have  been 
clearly  outlined,  but  it  remains  the  case  that  the  suitability  and  accuracy  of  the 
the  probes  will  be  a  matter  of  judgement  in  each  particular  application. 


5.  Conclusion 


A  design  for  a  thermocouple  probe  to  measure  solid  rocket  exhaust 
temperahires  has  been  produced.  A  method  of  correcting  for  conduction  losses 
in  the  probe  has  been  derived.  Both  the  probe  design  and  the  correction  were 
based  on  information  from  a  literature  survey  on  high  temperature 
measurement. 


PROBE  CALCULATIONS 


SUPERSONIC  STREAM  FLOW 
BOOST 


9/))/93 


Location 
Axial  distance 
Radial  distance 

PARTICULAR  FLOW  VARIABLES 


3 

0  096 


U1 

Ttofal 

Cp 

MW 


935 

2056 

1501 

28.45 


Gamma 

R 

1.2417602 

292.23199 

Trio/  stroam  tomporaturo 

Tstreom 

1765 

BEFORE  SHOCK 

Pressure 

PI 

Pol 

1.013E+05 

2.220E+05 

Temperoture 

T1  1765 

Tol  2056 

Mach  No 

Ml 

1.168 

AFTER  SHOCK 

Pressure 

P2 

Po2 

1  423E+05 

2.209E+05 

Temperature 
r2  1887 

To2  2056 

Mach  No 

M2 

U2 

0.860 

712 

IN  PROBE 

Pressure 

P3 

Po3 

2.086E-^05 

2.209E+05 

Temperature 

T3  2033 

To3  2056 

Moch  No 

^  A  -) 

U3 

0  3045 

262 

Expo 

Area  Ratio 

4  636 

2.000 

Density  Rho  0  351069 

(Adjust  M3  to  give  correct  areo  ratio) 

CALCULATIONS 

H  t  coef 

h 

2003.308 

D3 

Foci 

36.65015  D4 

355.2815 

2.278219 

LTFacl 

D5 

1.421126 

727.4386 


VELOCITY  CORRECTION  3 

CONDUCTION  CORRECTION  45 

TOTAL  ERROR  49 


» 


» 


CORMCTEO  TOTAL  TEMPERATURE 
CORRECTED  STREAM  TEMPERATURE 


2105 

1807 


VISCOSITY  AND  THERMAL  CONDUCIVITY  OF  GAS  IN  PROBE 


Gas 

Level 

Viscosity 

Thermal  conductivity 

Air 

0.65 

6  63E-05 

O.I09ft^4 

H20 

0.22 

7.4E-05 

0.258151 

C02 

0.13 

6.43E-06 

0.13557 

DISTIN 

UIN 

TOIN 

CPIN 

MWIN 

AIRIN 

H20IN 

C02IN 

SUMIN 

98 

935 

2056 

1501 

28.45 

21.0 

7.2 

4.2 

32.4 

FLOW  DEPENDENT  VARIABLES 


Viscosity  In  Probe 

Mu 

6.77E-05 

Thermal  conductivity  of  fluid 

Kf 

0.145606 

(Dofa  from  Rozrvevic) 

Thermal  conducivlty  of  wire 

Kw 

126.9675 

(Dota  from  Razrijevic) 

Biot  Number 

Bi 

0.04441 

Prondtl  Number 

Pr 

0.698314 

Nusselt  Number 

Nu 

6.879229 

Reynolds  Number 

Re 

677.7951 

I 


I 


PROBE  CONSTANTS 

Diameter  of  wire 

d 

0.0005 

Exposed  length  of  TC 

LI 

4.00E-C 

Total  length  of  TC 

L2 

5.40E-C 

Gas  constant 

Rbar 

8314 

Area  of  one  wire  Aw  1 .96E-07 

Exposed  Area  of  TC  Aex  1 .26E-05 

(2  wires) 


Figure  18:  Spreadsheet  output  of  calculation  of  conductio>i  correction 


» 
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